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Abstract 

We prove that the Schafer-Wayne short pulse equation (SPE) 
which describes the propagation of ultra-short optical pulses in non- 
linear media is integrable. First, we discover a Lax pair of the SPE 
which turns out to be of the Wadati-Konno-Ichikawa type. Second, 
we construct a chain of transformations which relates the SPE with 
the sine-Gordon equation. 



1 Introduction 

Recently, Schafer and Wayne [T] derived the short pulse equation (SPE) as a 
model equation, alternative to the nonlinear Schrodinger equation (NLSE), 
to approximate the evolution of very short optical pulses in nonlinear media. 
More recently, Chung, Jones, Schafer and Wayne [2j proved numerically that, 
as the pulse length shortens, the NLSE approximation becomes less and less 
accurate while the SPE provides a better and better approximation to the 
solution of Maxwell's equation. 

In the present paper, we study this interesting new nonlinear equation, 
the SPE, from the standpoint of its integrability. The NLSE is well known as 
one of the numerous nonlinear equations integrable by the inverse scattering 
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transform technique Therefore it is natural to ask whether the SPE, as 
an ultra-short pulse alternative to the NLSE, is integrable as well, or one has 
to rely mainly on numerical techniques when studying the SPE. We prove 
that the SPE is integrable. In Section El we find that the SPE possesses a 
Lax pair of the Wadati-Konno-Ichikawa (WKI) type jlj. Then, in Sectional 
we show how the SPE can be transformed into the well-known integrable 
sine-Gordon equation [3]. Section HI summarizes the results. 
We study the SPE in the form 

u xt = u+l(u 3 ) xx (1) 

which is related to the original form used in QjJ |2| by scale transformations 
of the dependent and independent variables. 



2 Zero-curvature representation 

Let us try to find a zero-curvature representation (ZCR) of the SPE (JTJ), 

D t X-D x T+[X,T} = 0, (2) 

which is the compatibility condition of the overdetermined linear system 

* x = XV, (3) 
% = Ttf , (4) 

where X and T are n x n matrix functions of u and its derivatives, ty(x,t) 
is an n-component column, D x and D t denote the total derivatives, and the 
square brackets denote the matrix commutator. 
Taking for simplicity 

X = u x A + B, T = T{u,u x ), (5) 

with constant matrices A and B, and replacing u x t by the right-hand side of 
flU), we find from © that 

T= ±u 2 u x A + S(u), (6) 

where the matrix function S(u) satisfies the equations 

dS/du = [A,S — \u 2 B] , uA + [B, S] = 0. (7) 



2 



Setting S(u) to be the simplest nontrivial polynomial, 



S = \u 2 B + uP + Q 



(8) 



with constant matrices P and Q, we have from (J7J) 
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[B,P] 



(9) 



and 



[P, [B, P}] = B 



[B,Q) = 0, 



[Q,[B,P]]=P 



(10) 



A nontrivial solution of the commutator equations (fTU|) can be easily 
found in traceless 2x2 matrices B, P and Q. Using this solution, together 
with ©, ©, © an d ©, and simplifying the expressions for X and T by 
X i — y GXG^ 1 and T i— > GTG^ 1 with a constant matrix G, we obtain the 
following: 



where A is an arbitrary nonzero constant. 

Consequently, the SPE (JTJ) possesses the ZCR (J2|), or the Lax pair (jSJ) 
and (J3J) , with the matrices X and T determined by and ()12j) . The linear 
problem (jSJ) with the matrix X (jllj) is a spectral problem of the WKI type 
j3] with respect to u x . Note, however, that the matrix T (|12|) is nonlocal 
with respect to u x and contains A -1 . The WKI- type ZCRs with matrices 
T containing negative powers of the spectral parameter were studied in j^J. 
Nevertheless, the nonlinear equation we study, the SPE (JTJ, did not appear 
in the literature prior to d 121; as far as we know. 

We have already pointed out that the form ((TJ) of the SPE is related to 
the original form used in |2] by scale transformations of variables. Note, 
however, that some complex- valued scalings of variables are required in order 
to change the sign at the nonlinear term in Therefore, if one does not 




(11) 



and 




(12) 
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accept complex- valued transformations (say, as nonphysical) but is interested 
in the form 




xx 



(13) 



of the SPE, one needs to use the ZCR ((21) with the matrices 





A 



\u x 
-A 



) 



) 



(14) 



for the 'minus' form (J13j) of the SPE. All this is completely analogous to the 
situation with the 'plus' and 'minus' forms of the NLSE [3]. 

3 Equivalence transformation 

Quite often an interesting newly-found nonlinear equation turns out to be 
equivalent to a well-studied old one through some chain of transformations. 
Let us see that the SPE (0) is not an exception in this respect. Of course, 
since the SPE possesses a WKI-type Lax pair, it is possible to use the known 
interrelation E| between spectral problems of the WKI type and the 
Ablowitz-Kaup-Newell-Segur (AKNS) type, in order to find which equa- 
tion possessing an AKNS-type Lax pair corresponds to the SPE. However, 
we follow a different way when trying to transform the SPE into some old 
equation: we use a generalized symmetry of the SPE as an indicator of re- 
quired transformations. This heuristic method is based on the facts that any 
generalized symmetry is formally represented by an evolution equation jS] 
and that the transformations of evolution equations are studied much better 
than transformations of other types of equations. 

It is easy to find directly that the SPE ((H) admits a third-order generalized 
symmetry, which is written in the evolutionary form as 



where the subscript s denotes the derivative with respect to the group pa- 
rameter s. This symmetry can also be obtained using the recursion operator 




(15) 
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R = Dl ■ u- l x D x ■ (ul + 1) 1/2 D- 1 ■ (ul + 1)" 3/2 u xx (16) 

easily derivable from the matrix X (|11|) by the cyclic basis technique [3 CD] . 
This operator generates the symmetry (|15|) from the trivial symmetry it s = 0, 
not from the translational symmetry u s = u x , because Ru x = 0. In what 
follows, however, we do not study the recursion operator (fTHjl. 

We consider the generalized symmetry ()15j) as a third-order evolution 
equation in u(x,t, s), and try to transform it into a simple known evolution 
equation, following the way used in [TO] HI] . 

First of all, we find the separant df /du xxx of the equation (fTHj) . where / 
denotes the equation's right-hand side, introduce the new dependent variable 
v(x,t, s) = (df/du xxx ) l ^ } and thus obtain the transformation 

v(x,t,s) = (u 2 x + l)- 1/2 (17) 
which relates the equation (JT5j) with the equation 

3 3v 2 v x v xx 3v 3 vl 

v s = v v xxx H — H 77 (18) 

1 _ v 2 ^ _ ^2)2 v ; 

possessing the separant v 3 . 

Then, we introduce the new independent and dependent variables, y and 
w(y, t, s), such that 

x = w(y,t,s), v(x,t,s)=w v (y,t,s). (19) 

This transformation (|T9"|) relates the equation (|T5j) with the constant-separant 
equation 

W S = Wyyy + ~ ^ + (f)(t, s) W y , (20) 

2 (1 - wf) 

where the arbitrary function <p(t, s) appeared as a 'constant' of the integra- 
tion over y. This inessential function </>(£, s) can be turned into zero by the 
redefinition of y, y i— > y + J (pit, s) ds, which does not change v(x, t, s). 

At last, trying to eliminate the second-order derivative terms from the 
transformed equation, we find that the new dependent variable 

z(y, t, s) = arccos(wj / ) (21) 
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transforms the equation (|2(J|) (where we have already set 4>(t, s) = 0) into the 
potential modified Korteweg-de Vries equation 



Z S — Zyyy + ^ Z y . (22) 

We have transformed the generalized symmetry (fTH|) of the SPE into 
the simple known evolution equation (j22j). The equation represents a 
generalized symmetry of the sine-Gordon equation z yt = sin z [5j, but not 
only of it: for example, the equations z yt = cos z, z yt = exp(±iz) and z yt = 
also admit the generalized symmetry (J22j) . In order to see which equation 
possessing the generalized symmetry (J22j) corresponds to the SPE, we apply 
the chain of transformations (|T7|). (fT§j) and (|7T|) directly to the SPE (JTJ), 
and actually obtain the sine-Gordon equation z yt = sin z. Note that, when 
making the transformation (jl9j) . we neglect the 'constant' of the integration 
over y, like we did this in the equation (|2U|) and for the same reason. 

Consequently, the SPE ((TJ) and the sine-Gordon equation 

z yt = sin ,2 (23) 
are equivalent to each other through the chain of transformations 

v(x,t) = (ul + l)- 1/2 ; (24) 
x = w(y,t), v(x,t) =w y (y,t); (25) 
z(y,t) = arccos(w y ). (26) 



4 Conclusion 

In this paper, we studied the integrability of the short pulse equation (SPE) 
which describes the propagation of ultra-short optical pulses in nonlinear 
media. We found a Lax pair of the SPE, of the Wadati-Konno-Ichikawa 
type. This makes possible to use the powerful inverse scattering transform 
technique for obtaining and analyzing solutions of the SPE. Also we found 
how to transform the SPE into the well-studied sine-Gordon equation. This 
makes possible to derive solutions and properties of the SPE from known 
solutions and properties of the sine-Gordon equation. We believe that these 
two possibilities to study the SPE analytically, discovered in this paper, will 
be valuable for soliton theory and nonlinear optics. 
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